Constraining Isocurvature Perturbations with CMB Polarisation 
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The role of cosmic microwave background (CMB) polarisation data in constraining the presence 
of primordial isocurvature modes is examined. While the MAP satellite mission will be unable to 
simultaneously constrain isocurvature modes and cosmological parameters, the PLANCK mission 
will be able to set strong limits on the presence of isocurvature modes if it makes a precise measure- 
ment of the CMB polarisation sky. We find that if we allow for the possible presence of isocurvature 
modes, the recently obtained BOOMERANG measurement of the curvature of the universe fails. 
However, a comparably sensitive polarisation measurement on the same angular scales will permit 
a determination of the curvature of the universe without the prior assumption of adiabaticity. 



Measurements of the cosmic microwave background 
anisotropy may soon allow us to rigorously determine 
the fundamental character of the primordial cosmological 
perturbations with a minimum of simplifying hypotheses. 
Recent data already impose formidable constraints 
on the parameter space of simple inflationary models, and 
with forthcoming satellite missions this situation is 
expected to improve substantially. 

This paper focuses on how to test the assumption of 
adiabaticity, namely that all components contributing to 
the density of the Universe are present in spatially uni- 
form ratios on hypersurfaces of equal cosmic temperature 
and initially share a common velocity field. The hypoth- 
esis of adiabaticity, put forth initially on the basis of sim- 
plicity, gained support when it was realized that this was 
the prediction of the simplest one- field inflationary iodels 
||. It should be noted, however, that multi- field infla- 
tionary models generically excite isocurvature modes as 
well @. 

In order to test the hypothesis of adiabaticity through 
observation it is necessary to study models where the pri- 
mordial perturbations are not solely adiabatic and thus 
to attempt to set bounds on the allowed admixtures of 
non-adiabatic modes. Non-adiabatic, or isocurvature, 
perturbations have already been studied in the litera- 
ture, but most work has focused on evaluating the via- 
bility of cosmological models in which the perturbations 
where entirely isocurvature in character, with no adia- 
batic component at all . Observational constraints on 
an uncorrelated admixture of adiabatic and cold dark 
matter isocurvature perturbations have also been con- 
sidered ||. 

Baryon isocurvature (BI) and cold dark matter isocur- 
vature (CDMI) models were studied some time ago 0. 
More recently it was realized that two additional isocur- 
vature modes are possible: a neutrino isocurvature den- 
sity (NID) mode and a neutrino isocurvature velocity 
(NIV) mode In the NID mode, the neutrino- 

photon ratio varies spatially. As modes enter the horizon, 
the photon-baryon fluid begins to oscillate acoustically 
whereas the neutrinos free stream. This differential be- 
havior perturbs the total energy density, leading to struc- 



ture formation via gravitational clustering. The neutrino 
velocity mode assumes a relative velocity between the 
photon and neutrino components but zero initial total 
momentum density. As with the NID mode, differential 
evolution spoils this cancellation after horizon crossing, 
again leading to structure formation. 

In a universe composed of just photons, baryons, neu- 
trinos, baryons, and a cold dark matter component, these 
four isocurvature modes and the adiabatic mode exhaust 
the possible modes nonsingular in the t — ► limit [||. 
The most general Gaussian primordial perturbation in 
such a cosmology is completely characterized by the ma- 
trix valued generalization of the power spectrum 

(A a (k) A b (k')) = P ab (k)-S 3 (k-k>) 

where the indices (a, b = 0,1,2,3,4) label the modes. 
When expectation values of observables quadratic in the 
linearized perturbations (such as the CMB Ci 's) are con- 
sidered, the assumption of Gaussianity is superfluous. A 
detailed discussion of how the CMB satellite missions will 
be able to constrain the presence of these modes is given 
in ref. Q. 

We briefly remark on the possible microphysical ori- 
gin of these modes. The neutrino isocurvature mode is 
rapidly damped before neutrino decoupling, and there- 
fore can be plausibly produced only by physics operat- 
ing after ~ 1 second. The neutrino isocurvature density 
mode is likewise damped by electroweak B + L violating 
anomalous processes (which would convert it mostly into 
a baryon isocurvature mode) operating at times earlier 
than 10~ 10 seconds. Since both times are well before pho- 
ton decoupling (~ 10 13 seconds), we nevertheless think it 
legitimate to describe the perturbations as 'primordial'. 

The temperature anisotropy spectra associated with 
the regular perturbation modes and their cross- 
correlations are shown in Fig. 1. An appropriate power 
law auto-correlation spectrum was assumed for each 
mode, so that the large scale CMB anisotropy is approxi- 
mately scale invariant. The cross-correlation power spec- 
tra were then taken to be proportional to the geometric 
mean of the two auto-correlation spectra. Clearly there is 
scope for significant generalisation of these assumptions. 
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1. CMB multipole spectra for the various modes, their cross-correlations, variations in the cosmological parameters. 
From top to bottom the rows show 1(1 + l)Ci/2n for the temperature, polarization, and temperature-polarization cross correla- 
tion, respectively, in fiK. The C; spectra for the various modes and their cross correlations are shown in the first two columns. 
The rightmost column shows the derivatives of the spectra with respect to the different cosmological parameters. The modes are 
indicated as follows: adiabatic (AD), neutrino isocurvature velocity (NIV), baryon isocurvature (BI), and neutrino isocurvature 
density (NID). A fiducial model with the parameter choices fi,;, = 0.06, J1a = 0.69, f2 C£ j m = 0.25, h — 0.65, r re i „ =0.1 and 
n s — 1 has been assumed. Because the CDM isocurvature mode produces a spectrum nearly identical to that of the BI mode, 
it is not considered separately. 



The adiabatic C'i temperature spectrum is charac- 
terised by a flat Sachs- Wolfe plateau at low I and a series 
of acoustic peaks, at I s» 220 + 310n, (n = 0, 1, . . .) in a 
flat Universe. Scale invariant baryon and CDM isocur- 
vature spectra produce little power at high I. The NID 
mode exhibits the phase shift characteristic of isocurva- 
ture density modes but the NIV mode produces a pattern 
of peaks more similar to that of the adiabatic mode. The 
NIV mode acquires a nearly cancelling phase shift be- 
cause the velocity is out of phase with the density. Note 
in particular how similar the NIV-adiabatic mode cross 
correlation C; is to the pure adiabatic mode Ci . The po- 
larisation and temperature-polarisation cross correlation 
power spectra associated with the isocurvature modes are 
indicated in the bottom two rows of Fig. 1. 

How feasible will it be to constrain or to detect isocur- 
vature modes using CMB measurements? The key ques- 
tion is whether one can distinguish the effect of the 
isocurvature modes from those of variations in the cosmo- 
logical parameters (see e.g. |^2|). The derivatives of the 
Ci power spectra about a fiducial ACDM model with re- 
spect to cosmological parameters are shown in the right- 
most column of Fig. 1. For small admixtures of isocurva- 
ture modes, the question is whether linear combinations 



of these spectra can be distinguished from those of the 
first two columns of Fig. 1. 

For small variations of cosmological parameters and 
small admixtures of isocurvature modes one can parame- 
terize the likelihood function as a multivariate Gaussian 
about a fiducial adiabatic model. Using projected esti- 
mates of the instrument noise for the MAP and PLANCK 
satellite missions, we have computed the estimated errors 
on the cosmological parameters and isocurvature auto- 
and cross-correlation amplitudes assuming the sky is ac- 
tually described by a simple adiabatic ACDM model. Ta- 
ble I shows that the MAP satellite will be unable to simul- 
taneously constrain isocurvature modes and measure the 
cosmological parameters. The PLANCK satellite will, 
but only if it measures the CMB sky polarisation as ac- 
curately as currently planned. 

We now discuss precisely how the polarisation mea- 
surement resolves the degeneracy between isocurvature 
modes and cosmological parameters. In Fig. 2 we illus- 
trate how polarization serves to remove the degeneracy 
corresponding to the eigenvector pointing in the most un- 
certain (i.e., flattest) direction of the relative likelihood 
when only temperature information is taken into account. 
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TABLE I. This table indicates the one sigma percentage errors on cosmological parameters and isocurvature mode amplitudes 
anticipated for the MAP and PLANCK satellite experiments. In the column headers, T denotes constraints inferred from 
temperature measurements alone, TP those from the complete temperature and polarisation measurements, and T+P those 
inferred if temperature and polarisation information is used separately without including the cross-correlation. 



The dotted curves are the contributions of the various 
components of the eigenvector. Summing these contri- 
butions, one finds very little net contribution to the tem- 
perature anisotropy (solid line) . This sum is shown mul- 
tiplied by ten for clarity. The lower two panels indi- 
cate the corresponding curves for the polarization and 
temperature-polarization cross-correlation, respectively. 
For these the delicate cancellation is broken, indicating 
that the polarization and temperature-polarization cross- 
correlation provide the information required to break the 
degeneracy. Fig. [| shows constraints on the four small- 
est temperature eigenvalues that will be provided by a 
PLANCK measurement of the polarisation and cross- 
correlation spectra. The area under the curve gives the 
Fisher information, or log likelihood in the Gaussian 
approximation, provided by the polarisation measure- 
ments. The degeneracy breaking by PLANCK polari- 
sation measurements occurs primarily at Z<100. Fig. 2 
shows that there is in fact considerable degeneracy break- 
ing at higher I, but this is not detectable with the instru- 
ment noise anticipated for PLANCK. 

To see why considerable information resides in the 
range Z<100, note that on these scales one directly ob- 
serves the primordial (super-horizon scale) polarisation. 
This is very different for the adiabatic and isocurvature 
modes (Fig. 4). In particular, for the NIV mode, since 
the large scale CMB anisotropy comes mainly from the 
Doppler effect rather than the Sachs- Wolfe effect, the po- 
larisation spectrum is enhanced by a factor of l~ 2 at low 
I relative to that for the adiabatic mode. The effect at 
very low I of varying T re i on , is well known fll2|] . 
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FIG. 2. Breaking the degeneracies with polariza- 
tion. The top panel indicates the delicate cancellation in 
the temperature power spectrum between the various com- 
ponents of the most uncertain principal direction. The lower 
panels show how this cancellation is broken in the polarization 
and temperature-polarization cross-correlation spectra. 
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FIG. 3. When only temperature information from 
PLANCK is taken into account, the uncertainties in the four 
most poorly measured principal directions are 239%, 60%, 
36%, and 23%. These numbers are the inverse square root of 
the corresponding eigenvalues of the Fisher matrix. When the 
polarization information anticipated from PLANCK is taken 
into account as well, these uncertainties are reduced to 11.1%, 
10.3%, 6.6%, and 4.6%, respectively. The plots (from top to 
bottom, respectively) indicate the contributions of the polari- 
sation information at each I to diagonal elements of the Fisher 
matrix in these directions. The cross correlation contribution 
is by definition the difference between the total and that ob- 
tained from polarisation and temperature information taken 
separately. 



We do not account for a calibration uncertainty and fix 
all cosmological parameters except fi^. We then allow 
arbitrary amounts of isocurvature and cross-correlation 
power and attempt to set limits simultaneously on these 
and on SQk- With the assumption of adiabaticity, the 
data yield a one sigma error on SQk of 2%. However, with 
isocurvature modes allowed, this error rises to 577%. Of 
course, the approximation that the likelihood is Gaus- 
sian breaks down at such a large level. Nevertheless, 
the conclusion that the error in SQ^ is of order unity 
is firm. If we assume that it will be possible to make 
a BOOMERANG-like polarisation map with accurate 
source subtraction and that the polarisation error is op- 
timal, (Tp = 2cr|n boom ' then we nn d that the additional 
polarisation information allows the constraint on to 
be reduced to 13%, a very significant improvement. This 
work therefore provides strong motivation for such a po- 
larisation measurement. 
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FIG. 4. Polarisation signal of the various modes on large 
angular scales (Z<100). Measurements on these scales are 
largely responsible for the degeneracy breaking which polari- 
sation measurements allow. 

Finally, one can ask to what extent the constraint 
on the curvature of the universe, derived using the re- 
cent BOOMERANG data, 5fl k < 0.12 at 2a @, is af- 
fected by the possible presence of isocurvature modes. 
Adopting the above flat fiducial model, we have com- 
puted the errors in Sflk allowing for sample variance and 
using the published values for the instrument noise M. 
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